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The nonequilibrium dynamic phase transition, in the kinetic Ising model in presence of an 
oscillating magnetic field, has been studied both by Monte Carlo simulation (in two dimension) 
and by solving the meanfield dynamical equation of motion for the average magnetization. The 

• temperature variations of hysteretic loss (loop area) and the dynamic correlation have been studied 
' near the transition point. The transition point has been identified as the minimum-correlation point. 

The hysteretic loss becomes maximum above the transition point. An analytical formulation has 
been developed to analyse the simulation results. A general relationship among hysteresis loop area, 
(— I I dynamic order parameter and dynamic correlation has also been developed. 

D ' PACS number(s): 05.50.+q 

; I. INTRODUCTION 

The dynamics of magnetization reversal in simple ferromagnetic systems has recently attracted considerable scientific 
' interest to study the nonequilibrium responses. In this regard, the dynamical responses of the Ising system in presence 
of an oscillating magnetic field have been studied extensively . The dynamical hysteretic response and the 
nonequilibrium dynamical phase transition [Q-^ are two main subjects of interest to study the dynamic responses of 
the kinetic Ising model in presence of an oscillating magnetic field. 
Q 1 Tome and Oliviera [Q first studied the dynamic transition by solving the mean field (MF) dynamic equation of 
^ ' motion (for the average magnetisation) of the kinetic Ising model in presence of a sinusoidally oscillating magnetic 
, field. By defining the order parameter as the time averaged magnetisation over a full cycle of the oscillating magnetic 
field, they showed that the order parameter vanishes depending upon the value of the temperature and the amplitude 
^ of the oscillating field. In the field amplitude and temperature plane they have drawn a phase boundary separating 
^\ dynamic ordered (nonzero value of order parameter) and disordered (order parameter vanishes) phase. They [J] 
have also observed and located a tricritical point (TCP), (separating the nature (discontinuous/continuous) of the 
r transition) on the phase boundary line. 

' Since, this transition exists even in the static (zero frequency) limit such a transition, observed from the solution 
of mean field dynamical equation, can not be dynamic in true sense. This is because, for the field amplitude less than 
the coercive field (at temperature less than the static ferro-para transition temperature), the response magnetisation 
varies periodically but asymmetrically even in the zero frequency limit; the system remains locked to one well of the 
■ free energy and cannot go to the other one, in the absence of noise or fluctuation. On the other hand, in presence of 
\ thermal fluctuations, in the static limit, the system can go from one well to another via the formation of nucleating 
I ■ droplets. Vanishingly small fleld is required to push the system from one to other well. Consequently, the dynamic 
■ O phase boundary collapses, in the presence of thermal fluctuations. 

C To study the true dynamic phase transition (which should disappear in the static limit) one has to consider the 

effect of thermal fluctuations. In this regard, Lo and Pelcovits flrst attempted to study the dynamic nature 
of this phase transition (incorporating the effect of fluctuation) in the kinetic Ising model by Monte Carlo (MC) 
simulation. However, they ||^ have not reported any precise phase boundary. Acharyya and Chakrabarti ^] studied 
the nonequilibrium dynamic phase transition in the kinetic Ising model in presence of oscillating magnetic field by 

• extensive MC simulation. They have drawn the phase boundary and located a tricritical point (as observed) on 
d ■ the boundary. It has been also observed jHj that this dynamic phase transition is associated with the breaking of the 

symmetry of the dynamic hysteresis (m — h) loop. In the dynamically disordered (value of order parameter vanishes) 
phase the corresponding hysteresis loop is symmetric, and loses its symmetry in the ordered phase (giving nonzero 
value of dynamic order parameter). They have 1^ also studied the temperature variation of the ac susceptibility 
components near the dynamic transition point. It has been observed ^] that the imaginary or lossy (real) part of the 
ac susceptibility gives a peak (dip) near the dynamic transition point (where the dynamic order parameter vanishes). It 
was concluded that this is a possible indication of the thermodynamic nature of this kind of nonequilibrium dynamical 
phase transition. 
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The statistical distribution of dynamic order parameter has been studied by Sides et al ||7|]. The nature of the 
distribution changes (from bimodal to unimodal) near the dynamic transition point. They have also observed 
that the fluctuation of the hysteresis loop area grows and becomes considerably large as one approach the dynamic 
transition point. 

The relaxation behaviour, of the dynamic order parameter, near the transition point (in the disordered phase), 
has been studied 1^ recently by MC simulation and solving meanfield dynamic equation. It has been observed 
that the relaxation is Debye type and the relaxation time diverges near the transition point. The 'specific heat' 
and the 'susceptibility' also diverge Q near the transition point in a similar manner with that of fluctuations of 
order parameter and fluctuation of energy respectively. These observations [|j (divergences of fluctuations) indirectly 
supports the earlier facts |^ where the distribution of the dynamic order parameter becomes wider and the fluctuation 
of hysteresis loop area becomes considerably large near the transition point. 

Recently the experimental evidence pO[ | of dynamic transition has been found. The dynamical symmetry breaking 
(associated to the dynamic transition) across the transition point of the hysteresis loop, has been observed, in highly 
anisotropic (Ising like) and ultrathin Co/Cu(001) ferromagnetic films by surface magneto-optic Kerr effect, as one 
pas ses through the transition point. The dynamical symmetry breaking in the hysteresis loops has also been observed 
PI in ultrathin Fe/W(110) film. However, the detailed natures of the dynamic transition and the phase boundary 
are not yet studied experimentally. 

In this communication, the dynamic phase transition has been studied in the kinetic Ising model in presence of a 
sinusoidally oscillating magnetic field by MC simulation and by solving the mean field dynamical equation of motion 
for the average magnetization. The temperature variations of the hysteresis loss (or loop area), the dynamic correlation 
and the phase lag are studied near the dynamic transition point. The paper has been organised as follows: in section 

II simple analytic forms are given for the loop area, dynamic correlation and dynamic order parameter. In section 

III a general relationship has been developed among the various dynamical quantities. In section IV the models are 
introduced and in section V the numerical results are given. The paper ends with a summary of the work in section 
VI. 



II. ANALYTIC FORMS OF THE LOOP AREA AND THE DYNAMIC CORRELATION NEAR THE 

TRANSITION POINT 

The form of the oscillating magnetic field is 

h{t) = hocosiut). (2.1) 

The dynamic order parameter is defined as 

Q = — J) m{t)dt, (2.2) 
27r / 

which is nothing but the time averaged magnetisation over a full cycle of the oscillating magnetic field. The hysteresis 
loop area is 



A 



■ j) mdh = h^uj <j) m{t)s\n{u)t)dt, (2.3) 



which corresponds the energy loss due to the hysteresis. The Dynamic correlation is defined as 

C =< m{t)h{t) > - < m{t) >< h{t) >, 

where < .. > denotes the time average over the full cycle of the oscillating magnetic field. Since < h{t) > = 0, one 
can write 



C = ^ ^ ra{t)h{t)dt ^ -^-^ ^ m{t) cos{ujt)dt. 



(2.4) 



The dynamic correlation has another physical interpretation. For the cooperatively interacting spin system, this is 
the negative of the time averaged spin-field interaction energy (per spin) (< Ef >— — f J^i '^i h{t)dt) over a 
complete cycle of the oscillating field. 

In the dynamically disordered (Q = 0) phase and near the transition point, the time series of the magnetisation 
{ra{t)) can be approximated as a square wave with a phase lag 5 with the applied sinusoidal magnetic field. 
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1 for <t <t/A + 5/uj 
m{t) = {-l for r/4 + (5/w < t < 3r/4 + ,5/w (2.5) 
1 for ?,t/A + 5/uj<t <2tt/uj, 

where t is the time period of the oscillating field and 5 is the phase lag between magnetisation m(t) and the magnetic 
field h{t) = ho cofi{ujt) . The value of the hysteresis loop area can easily be calculated as 

A = AhQsai{5). (2.6) 

This form of the loop area was also obtained 1^ from the assumption that it is approximately equal to four times the 
product of coercive field and remanent magnetization (here the remanent magnetisation equal to unity), where the 
coercive field is identified as /iosin(i5) (the change in field during the phase lag). Considering the same form of the 
magnetisation the dynamic correlation C can also be calculated exactly as 

C=^cos(5). (2.7) 

TT 



From the above forms of A and C it can be written as 



(4/io)2 (2/ioA)^ 



1. (2.8) 



The above relation tells that the loop area A and the dynamic correlation C is elliptically related to each other. It 
may be noted here, that the previously studied ac susceptibility components 1^ obey a circular relationship (x'^ + x"'^ 
= (too/^o)^), where mo is the amplitude of the magnetization. 

The ordered region (Q 7^ 0) can be approximated by considering the following form of the magnetization 

{ I for < t < r/4 + (5/cj 

m{t)^\ \-mr for t/4 + < t < 3t/4 + (2.9) 
[ 1 for 3r/4 + (5/tj < i < 27r/w. 

In the above simplified approximation, it was considered that the magnetisation can not jump to the other well, 
however the value of initial magnetisation is reduced by the amount m^. In the real situation it has been observed 
that this well is not fully square (as assumed above in the form o f m (t)), it has a cusp like (or parabolic) shape. For 



TO,. 2, the above functional form of m(t) will take the form of 2.5 and one can get the disordered (Q = 0) phase. 
Taking the above form of magnetisation the dynamic order parameter Q can be calculated as Q = (2 — mr)/2. It may 
be noted that, in this simplified approximation the dynamic order parameter Q is independent of phase lag S, which 
is not observed in the real situation (phase lag shows a peak at the transition point). However, this simple picture 
can anticipate the convex (looking from the origin) nature of the dynamic phase boundary. As the temperature 
increases to^ increases and it also increases as the field amplitude increases. In the simplest asumption, one can 
consider to^ is proportional to the product of ho and T. Demanding, rrir = 2 for the dynamic transition (Q = 0), one 
can readily obtain {hi^)dTd = constant. This equation tells that the dynamic phase boundary will be convex. The 
convex nature of the phase boundary remains invariant even if one assumes that is any increasing function of both 
T and (for example, power law; ~ T^h^, in this particular case the equation of the dynamic phase boundary 
becomes T^{hi^yji = constant, it is easy to see that this gives the convex shape of the dynamic phase boundary ). 
However, this very simple asumption can not describe the entire form of the phase boundary accurately, particularly 
near the end points ((/io)d = and Td = 0). 

III. GENERAL RELATION AMONG DYNAMIC ORDER PARAMETER, HYSTERESIS LOOP AREA 

AND THE DYNAMIC CORRELATION 

From the usual definitions (given in earlier section) of C and A, one can write 

1 /27rC A \ If,, 

I — 7— 1 = <i> m(t) exp[—iujt)at, 



V27r V ojho ujho J ^/27r 
where miuj) — -^^= ^ m{t) exp{^iujt)dt. So, 
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and 

A = —hoLuV^wlm {m{oj)). 

The general (complex) form of m{ui') will be 

m{u}') = \m{u)') \ ex.p{i(p) 

"^(^') = vfe + 4^)'^'^^P* [-*^^"' 

So, Q is related with A and C as follows 



27rCj 



Above equation gives the general relationship among Q, A and C. 

It has been observed that the steady response m(f), to a simisoidally oscillating magnetic field {h{t) = ha cos(a)t)), 
is periodic (with phase lag 6) and has the same periodicity (r = 27r/a;) of the field. So, one can write m{t) in a Fourier 
series as 

oo oc 

in{t) = Oo + a„ cos(na;t) + 6„ sin(na;t). (3.2) 

n=l n=l 

Prom the usual definitions of Q, ^ and C, it is easy to see that 
uo = Q, ai = 2C/ho and bi = A/{TTho). 
So, one can write 

2C A 

m{t) = Q+ — cos{ujt) + + —r- sm{ujt) + (3.3) 

ho irno 

Keeping only the first harmonic terms (ignoring higher harmonics) one can easily express the instantaneous magne- 
tization as 

m{t) =Q + mo cos{uJt - S) (3.4) 
where the amplitude of magnetization is mo = [(2C//io)^ + (^/(''"'lo))^]^''^ and the phase lag is 5 = tan~^(A/(27rC)). 

IV. THE MODEL AND THE SIMULATION SCHEME 
A. The Monte Carlo study 

The local field (at time t) at any site i, of a nc;arc;st neighbour ferromagnetic Ising model in the presence of a time 
varying external magnetic field h{t) with homogeneous and unit interaction energy can be written as 

hi{t) = J2<^j{t) + h{t) (4.1) 

j 

where ai{t) = ±1 and j runs over the nearest neighbour of site i. The local field (at site i) hi{t) has an external field 
part h{t), which is oscillating sinusoidally in time 

h{t) =hosm{2'Kft), (4.2) 

where ho and / are the amplitude and frequency of the oscillating field. 

According to heat-bath dynamics, the probability pi{t) for the spin ai{t) will be up at time t is given as 

^hiit)/KBT 

^ (,hiit)/KBT ^ Q-hi(t)/KBT ' (^•^) 
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where Kb is the Bohzmann constant which has been taken equal to unity for simpUcity. It may be noted here that 
the spin-spin interaction strength J has been taken equal to unity. The temperature T is measured in the unit of 
J /Kb- Field is measured in the unit of J. The spin ai{t) is oriented (at time t) as 

(T,(t + 1) = Sign[p,(i)-r,(t)] (4.4) 

where ri{t) are independent random fractions drawn from the uniform distribution between and 1. 

In the simulation, a square lattice (L x L) is considered under periodic boundary conditions. The initial condition 
is all spins are up (i.e., 'Ji{t = 0) = 1, for all i). The multispin coding technique is employed here to store 10 spins 
in a computer word consisting of 32 bits. 10 spins are updated simultaneously (or parallel) by a single command. 
All words (containing 10 spins) are updated sequentially and one full scan over the entire lattice consists of one 
Monte Carlo step per spin (MCSS). This is the unit of time in the simulation. The instantaneous magnetisation 
{m{t) — {^1 L'^)^iCri{t)) is calculated easily. Some transient loops were discarded to have a stable loop and all the 
dynamical quantities were calculated from the stable loop. 

This simulation is performed in a SUN workstation cluster and the computational speed recorded is 7.14 Million 
updates of spins per second. 



B. The Meanfield study 

The meanfield dynamical equation of Ising ferromagnet in the presence of time varying magnetic field is 

dm , / m{t) + h{t)\ ,^ . 

+ tanh(^^i^j, (4.5) 

where the external time varying field /i(i) has the previously described sinusoidal form. T is the temperature measured 
in the unit of zJ/Kb {z is coordination number and Kb is Boltzmann constant). This equation has been solved for 
'm{t) by fourth order Runge-Kutta method by taking the initial condition m{t = 0) = 1.0. The value of the time 
differential (dt) was taken lO^'^, so that the error is 0{dt^) ^ 10^^^. The frequency uj of the oscillating field is kept 
fixed {lj = 27r) throughout the study. Some transient loops were discarded and all the values of the response are 
calculated from a stable loop. 



V. RESULTS 



A. The Monte Carlo results 



In the MC simulation, a square lattice of linear size L = 1000 is considered. The frequency uj of the oscillating field 
has been kept fixed (w = 27r x 0.01) throughout the study. From the Monte Carlo simulation technique described above 
the m — h or hysteresis loops were obtained. Some (~ 600) initial transient loops were discarded to have the stable 
loop. From this one can easily estimate the length of the simulation. For the above choice of frequency, 100 MCSS 
are required to form a complete loop (or cycle), and 600 such loops were discarded. It has been checked carefully 
that the loop gets stabilised (within a reasonably useful errorbars) for this choice. The dynamic order parameter 
Q = § m{t)dt is readily calculated. The loop area A and the dynamic correlation C have been calculated from 
the usual definitions. The phase lag 5 (between field and magnetization) has been calculated by taking the difference 
between the positions of minimum of magnetization and the magnetic field All values of Q, A, 5 and C for a 
particular temperature were obtained by averaging over 10 different random samples to obtain the smooth variation. 
Fig. 1 demonstrates the dynamic transition (with dynamic symmetry breaking) and the related phemomena (e.g., 
temperature variations of A, S etc.) at a glance. For a fixed field amplitude ho — 0.7 the time variations of h{t) 
and m(t) are plotted for various temperatures in the pictures in left column and the corresponding m — h loops are 
shown in the right column. For very low temperature (topmost pictures of Fig. 1), since no spin flip occurs (within 
the time period) the magnetization m(t) remains constant (unity) and consequently the m — h loop is a straight line 
having zero loop area. The dynamic order parameter is unity. The concept of phase lag (between m{t) and h{t)) is 
not applicable here. Obviously the dynamic correlation is zero. After slight increase of temperature (pictures in the 
second row) some small number of spin flips occurs (within the time period). For some time, m(t) decreases from 
unity and again it becomes equal to unity. The phase lag is the frequency (lo) times the time difference between the 
positions of minimum of m,{t) and h{t). The m — h loop encloses a finite but small area giving Q less than unity. 
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The dynamic correlation starts to grow. As the temperature increases further the phase lag (5) and the loop area A 
increases (pictures in the third row) and the dynamic order parameter Q decreases. In all three cases, described so far, 
the asymmetric shapes of the m — h loops are observed due to asymmetric time variation of the response magnetisation 
m{t). The dynamic correlation decreases. As the temperature is very close to the dynamic transition temperature, 
(fourth row), where the time variation of the response magnetisation is almost symmetric giving maximum values 
of S. The m — h loop is symmetric and the dynamic order parameter Q is almost zero. The dynamic correlation 
C becomes negative and minimal. As one increases the temperaure further (last row), the phase lag decreases, and 
the loop area decreases. The dynamic correlation starts to grow further. It may be noted here that the conventional 
hysteresis ot m — h loop is observed in this region of temperature. As the temperature increases further the dynamic 
correlation grows, shows a maxima or peak and then decreases. The loop area monotonically decreases. 

The dynamical phase transition, via the dynamical symmetry breaking of the hysteresis loops, has been observed 
in highly anisotropic and ultrathin (2D Ising like) ferromagnetic films (Co/Cu(001) and Fe/W(110)) JlO| , pl| by using 
surface magneto optic Kerr effect study at room temperature. In the recent experimental study in ultrathin 
Fe/W(110), the dynamical symmetry breaking of the hysteresis loop was nicely depicted in Fig. 1 of Ref. 

The temperature variations of Q, S, C and A for two different values of field amplitudes ho are shown in Fig. 2. In 
both the cases, it has been observed that, near the dynamic transition point {Q — 0) the phase lag gives a peak and 
the dynamic correlation C gives a shallow dip. The dynamic correlation C gives a smeared peak much above (around 
T = 2.6) the static (ferro-para) transition point (Tc = 2.269..) (see Fig. 2). The hysteresis loop area A shows a peak 
above the dynamic transition point. 

It is possible to explain these observations from very simple analytical results described above (section II). The 
phase lag 6 becomes maximum near the dynamic transition point. So, according to the analytical formulation (for C 
and A) for a fixed value of the field amplitude as the temperature increases the loop area A{= AhQs'mS) starts to 
increases as the dynamic order parameter Q starts to decrease and above the dynamic transition point (complete spin 
reversal) the loop area will be maximum and after that A will start to decrease. Similarly, the dynamic correlation 
C will remain approximately equal to zero until a considerable amount of spin flip occurs and Q changes appreciably 
and then starts to increase. Above and near the transition point, where the phase lag S decreases as temperature 
increases and C = cos((5) should start to increase. Which has been observed indeed. However, near the transition 
point it gives a shallow dip, where the value of the dynamic correlation C in minimum and negative. The phase lag S 
should be less than or at most equal to tt/2. The field (h{t) = ho cos{ujt)) crosses zero first at the phase value tt/2 and 
it becomes minimum (maximum negative) at the value of phase (ujt) equal to tt. The response magnetisation, should 
change its sign (cross zero) within this period. This is true for uj ^ limit, however for finite but sufficiently high 
frequency, this will not happen. The phase difference more than 7r/2 would be observed yielding the unconventional 
shapes of m — ft, or hysteresis loops. In practice, it was observed that some asymmetric shape of the (m— h) loop gives 
the value of phase lag S slightly higher than 7r/2. In this region, cos((5 = 7r/2 + e) = — sin(e), which is negative and 
will show a shallow dip (cusp like shape) at the point where S is maximum. According, to the analytical prediction, 
the loop area A{— Ahosin{5)) should show maximum at the transition point. However, strictly speaking and in 
practice it has been observed that the loop area A becomes peaked above the transition temperature. Since the loop 
area is much more strongly dependent on the actual shape of the magnetisation (which is not a perfect square wave 
in the temperature range concerned here). As the field amplitude increases the transition points shift towards the 
lower temperature. The maximum of S also increases and consequently the dip of C becomes deeper and it remains 
negative over wider range of temperature (since S remains larger than 7r/2 over wider range). It may be noted that 
the dynamic correlation C becomes zero (in the disordered or Q = region) where the phase lag S = 7r/2 = 1.57080... 
The dynamic correlation C shows a smeared peak at quite higher temperature (above the Onsager value) , which was 
misinterpreted as a signature of the stochastic resonance. In the MF study (next section), it was shown that this 
is also present in the absence of fluctuations (or stochasticity) . 

A similar previous study showed that the ac susceptibility components would give peak (or dip) near the 
transition point. In that case, the susceptibility components were calculated from the phase lag 6. The phase lag S 
would show a peak at the transition point. As a consequence the susceptibility components would show peak (or dip) 
reflecting the behaviour of phase lag 6. However, in this case, the three measurements of phase lag S, loop area A and 
dynamic correlation C are completely independent, and indicate the transition point separately. 



B. The Meanfleld results 

By solving the above meanfleld equation the m—h or hysteresis loops were obtained. The dynamic order parameter 
Q = § m{t)dt is readily calculated. The loop area A and the dynamic correlation C have been calculated by using 
the above deflnitions. The phase lag (between held and magnetization) has been calculated by taking the difference 
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between the minima positions of magnetization and the magnetic field Q. Fig. 3 shows the temperature variations 
of Q, i5, G and A for two different values of field amplitudes /iq- In both the cases, it has been observed that, near 
the dynamic transition point {Q = 0) the phase lag gives a peak and the dynamic correlation C gives a shallow dip. 
The hysteretic loss A gives peak above the transition (dynamic) point. The dynamic correlation C gives a smeared 
peak much above (around T — 1.3) the static (ferro-para) transition point [Tc — 1.0) (for a closer view see Fig. 4). 

This high temperature peak of the dynamic correlation was misinterpreted as a signature of stochastic resonance 
[ p2t . This peak is indeed present in the case where the fluctuation is absent (MF case) . The appearance of this peak 
at higher temperature can be explained as follows: for much higher temperature the time variation of instantaneous 
magnetisation is no longer a square wave like and becomes almost si nuso idal with a phase lage 5. In a very simple 
view, it can be approximated as m{t) = moCOS^Lot — S) (from eqn 3.4; Q = at very high temperature). The 
dynamic correlation becomes C = m^is. cos((5). Where mg is the amplitude of the magnetisation which monotonically 
decreases as the temperature increases. The phase lag 6 also monotonically decreases in the higher temperature. 
Consequently cos((5) increases and toq decreases as temperature increases. So, one would obviously expect a peak at 
a finite temperature (high enough) where the competition, between fall of mo and rise of cos(5) with respect to the 
temperature T, becomes comparable. No stochasticity is involved in it! The loop area A also gives a peak above the 
transition point. Due to the similar reason given in the earlier section the dynamic correlation C gives a shallow dip 
near the transition point. 



VI. SUMMARY 



The dynamical response of the kinetic Ising model in presence of a sinusoidally oscillating magnetic field has been 
studied both by Monte Carlo simulation (in two dimension) and by solving the meanfield dynamical equation of 
motion for the average magnetization. 

A general relationship amo ng the hysteresis loop area A, dynamic order parameter Q and the dynamic correlation 



C has been developed (eqn. 3.1). The time series of the magnetization can be decomposed in a Fourior scries and 



the constant term is identifed as the dynamic order parameter Q, the amplitudes of first harmonic terms arc found 



to be related to the hysteretic loss (for sine term) and the dynamic correlation (for cosine term) (eqn. 3.3). 

The dynamic order parameter, the loop area and the dynamic correlation have been calculated separately (both 
from MC and MF study) and studied as a function of temperature. It was observed (in both cases) that the dynamic 
correlation shows shallow (negative) dip near the transition point. The dynamic transition point has been identified 
as the the minimum-correlation point. The hysteretic loss A becomes maximum above the dynamic transition point. 
In this sense, the dynamic transition point is not the maximum lossy point. It may be noted that the earlier study 
of the ac susceptibility suggests that the dynamic transition point would be the maximum-lossy point, since the 
imaginary part (or lossy- part) of the ac or complex susceptibility also shows a peak near the dynamic transition 
point. However, there is a remarkable distinction from the present study. In the earlier study the phase lag was 
calculated from the simulations and the ac susceptibility components were calculated from the phase lag. So, it is 
expected that the temperature variations of the phase lag will be refiected directly in the temperature variations of 
ac susceptibility components. But, in the present study the measurements of phase lag, dynamic correlation and the 
loop area are completely independent. This behaviour of the dynamic correlation is explained from a simple square 
wave like time variation of the instantaneous response magnetisation. The oversimplified asumption is incapable of 
explaining the peak position (above the transition point) of the hysteretic loss A. However, this simple picture can 
qualitatively describe the nonmonotonic temperature variations of A and C. 

The high temperature (above the static critical point Tc) peak of the dynamic correlation was misinterpreted [T^ ] 
as a signature of the stochastic resonance. This was also discussed and an analytical form of the dynamic correlation 
was proposed to show that the high temperature peak of the dynamic correlation is present even in the absence of 
fluctuations (or stochasticity). 

Along with the dynamic correlation, the dynamic transition can be identified by various thermodynamic quantities 
like ac susceptibility Q| , relaxation time , specific heat j|] , susceptibility and the fluctuations of dynamic order 
parameter and energy All these quantities indicate the thermodynamic natures of this kind of noncquilibrium 
dynamic phase transition by showing peak, dip or divergence near the transition point. 

The related phenomena of this kind of nonequilibrium dynamic phase transition in the kinetic Ising model are 
mostly based on observations and not yet analysed by using rigorous theoretical foundations of equilibrium statistical 
mechanics available so far. The experimental evidences ||l^,|ll| are still in the primitive stage. Experimentally, only 
the dynamic symmetry breaking of the hysteresis loops is observed near the transition point. However, the details 
study of the nature of the transition, the phase boundary and the associated phenomena (described above) has not 
yet done experimentally. 
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Figure Captions 



Fig.l. A pictorial demonstration of dynamic transition and associated phemomena. The figures in the left column 
represents the time variation of h{t) and m{t) for different temperatures and the corresponding m — h loop are shown 
in right column. Temperature increases from top to bottom. Monte Carlo results for L = 1000, w = 27r x 0.01 and 
ho = 0.7. 

Fig. 2. The Monte Carlo results for the temperature variations of Q, S, C and A for two different values of field 
amplitudes. Q (solid lines, (I) for ho = 0.9 and (II) for ho = 0.7), S (x for ho = 0.9 and O for ho = 0.7), C (A for ho 
= 0.9 and + for ho = 0.7) and A (★ for ho = 0.9 and □ for ho = 0.7). 

Fig. 3. The mean-field results for the temperature variations of Q, d, C and A for two different values of field 
amplitudes. Q (solid Uncs,(I) for ho = 0.3 and (II) for ho = 0.2), 5 (x for ho = 0.3 and O for ho = 0.2), C (A for ho 
= 0.3 and + for ho = 0.2) and A (* for ho = 0.3 and □ for ho = 0.2). 

Fig.4. The closer view of Fig. 3. for the dynamic correlation plotted against the temperature (for fixed field amplitude 
ho = 0.2). 



8 



Fig. 2 (PRE, Acharyya) 




Fig. 3 (PRE, Acharyya) 




Fig. 4 (PRE, Acharyya) 




